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ABSTRACT 

Dark matter may consist of weakly interacting elementary particles or of macroscopic 
compact objects. We show that the statistics of the gravitational lensing of high redshift 
supernovae strongly discriminate between these two classes of dark matter candidates. We 
develop a method of calculating the magnification distribution of supernovae, which can be 
interpreted in terms of the properties of the lensing objects. With simulated data we show 
that £ 50 well measured type la supernovae ( Am ~ 0.2 mag ) at redshifts ~ 1 can clearly 
distinguish macroscopic from microscopic dark matter if fl a }t 0.2 and all dark matter is in one 
form or the other. 



Subject headings: cosmology: theory, dark matter, gravitational lensing 



1. Introduction 

The nature of dark matter (DM) poses one of the most outstanding problems in astrophysics. There 
are essentially two alternative hypotheses. The dark matter may be microscopic, consisting of weakly 
interacting particles (WIMPs) such as SUSY neutralinos or axions, or else be macroscopic, compact 
objects such as primordial black holes (PBHs), brown dwarfs or old white dwarfs (MACHOs). Big Bang 
nucleosynthesis (BBN) puts a bound on the density in baryonic matter of il^h 2 £ 0.02 (or ^ 0.03 if one 
allows for inhomogeneous BBN), but the density of PBHs is not well constrained. It is possible that some 
hitherto unknown mechanism allows for DM that is dominated by macroscopic objects. For these reasons 
direct observational constraints on macroscopic DM of any density are very important. 

We propose a simple test for distinguishing macroscopic from microscopic dark matter. In this letter 
we consider only the opposing hypotheses that one or the other dominates. If the DM is microscopic, the 
clustered component, in halos, lenses high redshift supernovae (SNe). If the DM is macroscopic, most light 
beams do not intersect any matter - no Ricci focusing - and the SN brightness distribution is skewed to an 
extent that can be quantitatively distinguished from halo lensing. 



2. Properties of the Magnification Probability Distribution Function 

In this paper we consider the lensing of distant supernovae by discrete "lenses" . A lens is the smallest 
unit of mass that acts coherently for the purpose of lensing. This could be a galaxy halo or it could be a 
high mass dark matter candidate such as a PBH. 

We make the distinction between macroscopic and microscopic DM more quantitative by considering 
two mass scales. The first is defined by the requirement that the projected density be smooth on the scale 
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of the angular size of the source. This gives a maximum mass of 



GMU 1 ", (A) (),/,->/ 



(1) 



where A s is the physical size of the source and / is a geometric factor of order unity. If the unit of DM is 
smaller then this it is microscopic DM. Another, larger mass scale is defined by the requirement that the 
angular size of the source be small compared to the Einstein ring radius so that it can be considered a true 
point source if 
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where D s is the angular size distance to the source. If a lens is much below this mass the high magnification 
tail of the distribution function will be suppressed and the rare high magnification events will become 



time-dependent (Schneider & Wagoner 1987). The measured velocity of the expanding photosphere of a 
type la SN is around 1.0 - 1.4 x 10 4 km s" 1 flPatat et. al 1996| )which means X s ~ 40 - 57 At AU/week. 
The SN reaches maximum light in approximately one week. 

The background cosmology will be taken to be the standard Friedman-Lemaitre-Robertson- Walker 



i(t) 2 {dx 2 + D(x) 2 d£l) where the comoving angular size distance 
{i?sinh(x/i?), X) ^ sm (x/^)} (R = \H V^ — £1 — ^a| _1 ) for the open, flat and closed global 



(FLRW) with the metric ds 2 = dt 2 
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geometries respectively. Another relevance angular size distance is the Dyer-Roeder or empty-beam 



distance, D(x) ( |Dyer fc Roeder 1974| , |Kantowski (1998)| , note difference in notation) which is angular size 
distance for a beam that passes through empty space and experiences no shear. 



2.1. Magnification by a single lens 

Consider a single lens at a fixed coordinate distance from Earth. The path of the light is described by 
either of two lensing equations: 

r x = y-a(y,D l ,D s ) (3) 
rx = [l-/c i (x.)]y-a(y J A,-D.) (4) 

where is the positions of the lens relative to the undeflected line of sight to the source, y is the position 
of its image in the same plane and a is the deflection angle times the angular size distance. In equation (^) 
a negative background convergence, Kb, is included to account for the lack of background mass density that 
is assumed when D is used instead of D. Two magnifications, fx and ll, can be defined using equations (||) 
and (§) respectively. The requirement that the two lensing equations agree on the true size of an object 
results in the relation D{\) — [1 — K b{x)]D{x)- The explicit form of Kb(z) can be found by comparing the 



standard FLRW expression for D{\) with the solutions for D(x) found in Kantowski (1998) 



The probability that the lens is located between r±_ and r±_ + dr± is p(rj_)drj_ oc r±drj_. If the 
lens is spherically symmetric and the magnification is a monotonic function of r± the expression for the 
magnification can be inverted (at least numerically) to get r±(fi, D, D s ). Then probability of a lens causing 
the magnification 1 + 8ll can be found by changing variables. Lenses might also have properties such as 
mass, scale length, etc. which need to be averaged. 



For the case of a point mass lens the total magnification of both images is given by \x = f 2 + 2/ r\Jr 2 + 4; 
f = r±/R e (m, D, D s ). The Einstein radius of the lens is given by R 2 = AGmDiDi s / D s . The single lens 
distribution function is then 

p(Sfi)dSLi cx [(1 + Sfi) 2 - l] ~ 3/2 dSfi. (5) 
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The probability in is not normalizable; it diverges at small 8jl. This can be handled by introducing a 
cutoff in either Sfl space or in r±. The nature of this cutoff is not important as long as it is at sufficiently 
small 5fi or large r±. This will be clear when the total magnification distribution due to multiple lenses is 
considered. 

If the dark matter consists of microscopic particles clumped into halos, the entire halo will act as a 
single lens. In this case the Ricci focusing contribution to the magnification strongly dominates over shear 
distortions produced by mass outside of the beam ( Holz fc Wald (1998)| , Premadi, Martel fc Matzner 1991- ) 
and is then a function of only the local dimensionless surface density, «(y). Furthermore the lensing of the 
great majority of SNe will be quite weak which allows us to confidently make the linear approximation: 
5/i = 2[/«(y, Di, D s ) + Kf,]. This assumption has been well justified by many authors and will be confirmed 



by results in § 2.2 



For the purposes of this paper it will suffice to use a simple model for the surface density of halos. We 
use models with surface densities given by 
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This model behaves like a singular isothermal sphere out to y_i 



where it is smoothly cut off. 



In the following calculations, each halo is assumed to harbor a galaxy. At all redshifts a Schechter 
luminosity function fit to local galaxies is assumed with a = —1.07 and <p* = 0.01(1 + z) 3 h 3 Mpc -3 . The 
luminosities are then related to the circular velocity, V c , by the local Tully-Fisher relation, V c = V^L/L*) 0,22 
where V* = 200 km s _1 . The scale lengths are related to the luminosity through r c = r*(L/ X*) 1 / 2 with 
r* = 220 kpc. The precise values used for these parameters do not have a significant effect on the results of 
this paper. 



2.2. Total Magnification 



The total magnification of a source includes contributions from all the lenses surrounding the light 
path. To find the true path connecting a source to us, the lensing equation must be solved with multiple 



deflections (see Schneider, Ehlers & Falco (1992)). The magnifications due to different lens planes are in 
general nonlinearly coupled. However, if the deflections due to no more than one of the lenses is very 
weak the coupling between lenses can be ignored and their magnifications, (5/x or 5p,, will add linearly. 
This is a good approximation for the vast majority of light paths in realistic models. The validity of this 



assumption will be justified by the results and is further investigated in Metcalf (1999) . Furthermore 
numerical simulations and analytic arguments show that for both kinds of DM it is a good approximation 



to take the lenses to be uncorrelated in space (see Holz & Wald (1998), Metcalf 1998b). If in addition we 
take the lenses internal properties to be uncorrelated, the probability that the total magnification, 6p, s , of a 
point source is between Sfl s and 5jl s + dSfl s is 

N 

P(6p, s )ddp, s = ddp, s / Y[ [dSjii p(Sp,i)} 
J i=i 

N 
i=l 

where the Sfli is the contribution of the zth lens. 



The magnification <5/i s is defined as the deviation of the luminosity from its mean value. As a 
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Fig. 1. — Histograms representing the total magnification probability distribution for macroscopic DM and 
microscopic DM clumped into halos. The means of all the distributions are zero. For the macroscopic DM 
case all the matter in the universe is in the lenses. The shape of the distribution for DM halos is dependent 
on both the cosmology and the specific halo model assumed. This is a representative sample. 



result, the mean of the distribution P(S/i s ) must vanish.^] This combined with the requirement that both 
magnifications agree on the true size of a source results in the expression 1 — /t&(x) = (A(x)) 1 ^ 2 - I n this way 
the value of Kb(x) can be found by calculating the mean of (Q) numerically and a consistency check of the 
calculations can be made by comparing the results with the explicit values for D(x) and D(x)- These values 
agree to a few percent which is consistent with the uncertainty introduced by the discrete nature of the 
numerical calculation in the power law tail of the distribution. The minimum magnification, 6/j, m in, in the 
single lens distribution is set low enough that the resulting total distribution is independent of the cutoff. 

Figure [l] shows some examples of histograms made by producing random values Sjli drawn from 
the single lens distributions and then adding them to get the total magnification. The macroscopic DM 
distributions shown in figure [l] are independent of the lens mass and peak well below their mean and 
near the empty-beam solutions (corresponding to 8^, = —0.21, —0.12 and —0.084) because in these cases 
most lines of sight do not come very close to any lens. The probability that there are two lenses which 
individually give magnifications greater than Sfj, becomes appreciable only below the peak. This supports 
our approximation that whenever the lensing is strong it is dominated by one lens and the coupling between 
lenses is small at this redshift. In addition we have compared our results with the numerical simulations of 



Holz fc Wald (1998) and found excellent agreement. 



3. Distinguishing Dark Matter Candidates 

The apparent luminosity of a SN, l b, after lensing can be expressed in terms of either of the two 
magnifications, l a b = fil = jj,l/(£i). We wish to infer via the measured luminosities of a set of SNe, each 
located at a different redshift, from which distribution the magnifications were drawn and in this way 



2 The actual mean angular size distance should be slightly larger than the FLRW value because galaxies obscure some sources. 
Galaxies are presumably correlated with high density regions through which the magnification would be above average were 
they transparent. 
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surmise which DM candidate is most likely. To establish some insight into the magnitude of this effect, the 
differences in magnitudes between the average and the empty-beam solutions at z = 1 are: —0.25 mag for 
Q = 1, —0.14 mag for flat £! = 0.3 and —0.10 mag for open f2 = 0.3. 

Let us denote the probability of getting a data set, {Sfi}, given a model - either microscopic or 
macroscopic DM - as P({Sfi}\model) = J7- P(Sfii\model)d5^ii where the product is over the observed 
SNe. The model here includes sources of noise. This probability can be calculated numerically from the 



probability distributions discussed in § |2.2| . Because of Bayes' theorem we know that the ratio of these two 
probabilities is equal to the relative likelihood of the models being correct, the odds, given a data set. It is 
convenient to modify the odds into the statistic 



M p = —— In 



J dfl dri\p(il , £Ia)P ({£/•*} | macroDM, noise) 
J <if2 o <if2Ap(0 o , $}a)P halos, noise) 



(8) 



where p(£l ,Q\) is the prior distribution for the cosmological model based on independent information or 
prejudice. The measured M. v is expected to be large if DM is macroscopic and smaller if DM is microscopic 
or nonexistent. 

For the left hand plot in figure || five thousand simulated data sets were created, M p is calculated for 
each of them and their cumulative distributions plotted. The noise included in the simulation originates 
from both the intrinsic distribution of SN luminosities, presently corrected to ~ 0.12 mag, and the 
observational noise, presently an additional ~ 0.08 mag. For the left hand plot the noise is taken to be 
Gaussian-distributed in magnitudes with a standard deviation of 0.16 mag except for the dot-dashed curves 
which have Am = 0.2 mag. The cosmology is fixed in this plot, ie p(Q ,Q\) is a (5-function. M p can be 
calculated for a given data set and compared to this plot to determine its significance. It can be seen here 
that for 51 SNe (solid curve) at z = 1 the two distributions overlap at the 4% level, ie 96% of the time one 
of the DM candidates can be ruled out at better than the 96% confidence level. One of the advantages 
of M p is that it is close to Gaussian distributed with a mean that is independent of the number of SNe 
observed. In this way, once the cosmology and noise model is fixed, the value of M. p is a direct prediction 
of the kind of DM. 

The middle plot in figure illustrates the importance of some possible systematic uncertainties that 
arise from not knowing precisely the distribution of the noise. The solid curves are the same as in the 
left hand plot. The dotted curve is the extreme case where the noise is actually Gaussian distributed in 
magnification (there is a low magnitude tail), but M. v is calculated under the same assumptions as in the 
left hand plot. The dashed line in this plot is the case where the standard deviation is overestimated to be 
Am = 0.2 mag but is really Am = 0.16 mag. These errors in the noise model do not destroy the efficacy 
of the test, but they could be important if a long tail exists in the intrinsic distribution of luminosities and 
they become more important for smaller f2 D and J7a- 

The right hand plot in figure ^ addresses the question of differentiating between DM candidates 
without assuming specific values for the cosmological parameters, thereby making the conclusion cosmology 
independent. Here the prior is taken to be p(£l , ^a) = S(l — f2 Q — Q\) within a range in Q Q (Afl a = 0, 0.1 
and 0.2) centered on 0.3 and zero otherwise. The simulated data is the same here as for the solid curves in 
the two left hand plots. However, the integrations in (|^) would be prohibitively time consuming if the entire 
magnification distribution function were calculated for each trial cosmology. To simplify the calculation 
without loosing much of the test's effectiveness we use approximate, analytic test distribution functions. 
For the macroscopic DM case we use (|^) with the low magnification cutoff which insures that it gives the 
correct mean. Comparison of this approximation with the full multi-lens distribution shows that it is a good 
approximation especially for low o . For the microscopic DM/halo case we approximate the distribution as 



a Gaussian with an appropriate width (see Metcalf 1998a). This plot shows that not only is this simplified 



calculational technique adequate, but that one does not need to assume a precise cosmological model 
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Fig. 2. — Differentiating dark matter candidates: The cumulative distributions of the statistic A4 P . The 
cases where the true DM is macroscopic rise toward the right and microscopic DM cases rise toward the left. 
All the observed SNe are taken to be at z = 1 and in all cases S1 G + £1 a = 1 • On the left the different curves 
are for different numbers of observed SNe as marked with Am = 0.16 mag in all cases except the dot-dashed 
curves. The halos have total density f2^ = 0.27. The sold curves in the center panel are the same as in the 
left hand panel and the other curves are described in the text. In the right hand panel there are 75 SNe 
observed. The widths of the prior distributions are marked. 



to differentiate between DM candidates. Increasing the width of the prior beyond AS7 = 0.2 does not 
make much difference. The reason for this is that if the assumed cosmological parameters are significantly 
different than the true ones the distribution will be shifted to an extent that it is no longer consistent with 
the data. This shift would be confused for a lensing effect if the two kinds of distributions, illustrated in 
figure [I], were translations of each other, but they are not, even after noise is added. For the two DM cases, 
the modes of the magnification distributions follow different m—z relations, but their means are the same. 
For a fixed redshift, it is the distribution of luminosities about the mean that distinguishes the two cases. 

For open models = 0) it is more difficult to differentiate the DM candidates, but even in this case 
with 51 SNe at z = 1 and fl = 0.3 we expect to get better than 90% confidence at least 90% of the time. If 
S7 = 0.1 BBN constraints just barely allow for all DM being made of baryonic objects. In this case similar 
bounds to those shown in figure | for 51 SNe can be achieved with 200 SNe. However the means of the M v 
distributions are closer together in this case, making the test more susceptible to systematic errors in the 
assumed noise model. 

The power of lensing to differentiate DM candidates comes mostly from its ability to identify 
macroscopic DM. A positive detection of the lensing by microscopic DM halos will take more SNe, as will 
constraining the precise fraction DM in macroscopic form, unless correlations between SN luminosities and 



foreground galaxies are utilized (Metcalf 1998b, Metcalf (1999)) 



4. Discussion 



One concern in implementing the test described here is the possibility that type la SNe and/or 
their galactic environments evolve with redshift. This is also a major concern in cosmological parameter 
estimation from SNe. So far there is no indication that the colors or spectra systematically change with 
redshift ( Pcrlmuttcr et al. 1997 , Riess et al. 1998 ). Since the evolution of the magnification distribution is 
determined by cosmology it is in principle possible to make an independent test for systematic evolution in 
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the distribution of SN luminosities. 

Microscopic DM does not need to be clustered for this test to work. The clustering is added to make 
the calculations realistic. Clustering the microscopic DM to a greater or lesser extent would affect our 
results quantitatively, but the test would still be viable in more extreme cases. We conclude that if the 
assumptions we have made about the noise levels in future SN observations remain reasonable, on the order 
of 50 — 100 SNe at z ~ 1 should suffice to determine a fundamental question: whether the major constituent 
of extragalactic DM is microscopic particles or macroscopic objects. 



We would like to thank D. Holz for providing the results of his simulations for the purposes of 
comparison. 
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